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ABSTRACT Distinct lineages of avian influenza viruses (AIVs) are harbored by spatiaUy segregated birds, yet significant surveil- 
lance gaps exist around the globe. Virtually nothing is known from the Antarctic. Using virus culture, molecular analysis, full 
genome sequencing, and serology of samples from Adelie penguins in Antarctica, we confirmed infection by HI 1N2 subtype 
AIVs. Their genetic segments were distinct from aU known contemporary influenza viruses, including South American AIVs, 
suggesting spatial separation from other lineages. Only in the matrix and polymerase acidic gene phylogenies did the Antarctic 
sequences form a sister relationship to South American AIVs, whereas distant phylogenetic relationships were evident in all 
other gene segments. Interestingly, their neuraminidase genes formed a distant relationship to aU avian and human influenza 
lineages, and the polymerase basic 1 and polymerase acidic formed a sister relationship to the equine H3N8 influenza virus lin- 
eage that emerged during 1963 and whose avian origins were previously unknown. We also estimated that each gene segment 
had diverged for 49 to 80 years from its most closely related sequences, highlighting a significant gap in our AIV knowledge in 
the region. We also show that the receptor binding properties of the H11N2 viruses are predominantly avian and that they were 
unable to replicate efficiently in experimentaUy inoculated ferrets, suggesting their continuous evolution in avian hosts. These 
findings add substantiaUy to our understanding of both the ecology and the intra- and intercontinental movement of Antarctic 
AIVs and highlight the potential risk of an incursion of highly pathogenic AIVs into this fragile environment. 

IMPORTANCE Avian influenza viruses (AIVs) are typicaUy maintained and spread by migratory birds, resulting in the existence of 
distinctly different viruses around the world. However, AIVs have not previously been detected in Antarctica. In this study, we 
characterized HI 1N2 viruses sampled from Adelie penguins from two geographicaUy different sites in Antarctica and show that 
the segmented AIV genome diverged between 49 and 80 years ago from other AIVs, with several genes showing similarity and 
shared ancestry with H3N8 equine influenza viruses. This study provides the first insight into the ecology of AIVs in Antarctica 
and highlights the potential risk of an introduction of highly pathogenic AIVs into the continent. 
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WUd aquatic birds such as dabbling ducks, gulls, and other 
shorebirds are considered the natural reservoir for avian 
influenza viruses (AIVs). The ecology and migratory patterns of 
these birds therefore have a direct effect on the global distribution 
and diversity of AIVs (1). The ecology of AIVs in wild birds has 
been well studied in many regions of the Northern Hemisphere, 
but considerably less is known in the Southern Hemisphere. Of 
the publically available AIV sequences (n = 19,784), only 5.7%, 
1%, and 0.1% are from Africa, Oceania, and South America, re- 
spectively, and none have previously been described from Antarc- 
tica (GISAID, http://www.gisaid.org) . Wild migratory birds play a 
key role in the spread of AIVs on a local, regional, and interconti- 
nental scale via broadly established flyways (1) (see Fig. SI in the 
supplemental material). As a result, spatiaUy segregated birds of 



the American and Eurasian landmasses harbor distinct lineages of 
AIVs (1). Traditionally, AIVs have been broadly separated into 
North American or Eurasian lineages, although this appears to be 
an oversimplification and less applicable to AIVs from the South- 
ern Hemisphere (2). Recent genetic analysis of AIVs in Australia, 
located at the southern end of the East Asian-Australian flyway, 
showed some divergence from Eurasian strains, demonstrating 
that, although new strains have been introduced frequently from 
Europe and Asia, there is a high degree of endemic evolution that 
has produced a group of viruses that are genetically distinct from 
those circulating in Asia (2-4). Similar evolutionary divergence 
has also been observed in some South American AIVs, where in- 
ternal gene segments form unique "South American-like" clades, 
although the hemagglutinin (HA) and/or neuraminidase (NA) 
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TABLE 1 Prevalence of influenza virus and influenza antibodies in Adelie penguins 



Cloacal/traclieal swabs Serum samples 







No. of 


No. (%) of samples 


No. of samples 


No. of 


No. (%) of samples 


No. (%) of samples with 






penguins 


with AlV RNA 


isolated in 


penguins 


with lAV antibodies 


A(H11N2) antibodies 


Location 


Age 


sampled 


detected" 


eggs* 


sampled 


detected^ 


detected'' 


Admiralty Bay 


Adult 


110 


0 


0 


96 


16 (16.7) 


0 




Chick 


40 


1 (2.5) 


0 


35 


0 


0 


Rada Covadonga 


Adult 


111 


6 (5.4) 


3« 


99 


27 (27.3) 


2(1.8) 




Chick 


40 


1 (2.5) 


If 


40 


0 


0 


Total 




301 


8 (2.7) 


4 


270 


43 (15.9) 


2 (0.7) 



" No. of samples where the HA, NS, or M gene sequence was obtained and confirmed to be influenza A virus by BLAST analysis. 
^ Egg isolation was attempted only for those samples with a C^- value of <38 (n = 14). 

Determined using an anti-NP ELISA assay. 
^ Determined using an HI assay with the HI 1N2 antigen A/peng/ Ant/270/13. 

"Viruses were designated A/ Adelie Penguin/ Antarctica/226/2013, A/Adehe Penguin/An tarctica/184/2013, A/Adelie Penguin/ Antarctica/178/2013. 
/Virus was designated A/Adelie Penguin/Antarctica/270/2013. 



gene segments are often "North American-like" (5-8). It is there- 
fore possible that evolutionary divergence of AIVs may increase 
with further advancement south and that spatially segregated re- 
gions and continents such as Southern Patagonia and Antarctica 
may act as evolutionary "sinks" harboring highly diverged AIVs. 

The major migratory flyways are not thought to extend to Ant- 
arctica, although each spring, over 100 million birds breed around 
the rocky Antarctic coastline and offshore islands (9). These in- 
clude birds such as the Arctic tern (Sterna paradisaea) (10) and 
south polar skua (Stercorarius maccormicki) (11) that conduct 
transhemispheric migrations to Antarctica and thereby have the 
potential to facilitate movement of avian influenza viruses into 
and out of Antarctica. A small number of earlier AIV studies in 
Antarctica have detected influenza A antibodies in penguins and 
other birds (12-16), but none have detected or isolated any influ- 
enza viruses. The detection and subsequent analysis of AIVs from 



Antarctica would provide insight into the frequency and likely 
routes of AIV movement into and within the continent and help 
assess the potential risk of an incursion of highly pathogenic AIVs 
from neighboring regions. 

RESULTS 

Identification of influenza virus in penguins. We took combined 
cloacal/tracheal swabs from 301 Adelie penguins (Pygoscelis ade- 
liae) (adults and chicks), and blood from 270 of those, from two 
locations on the Antarctic Peninsula during January and February 
2013 (Table 1 and Fig. 1). Eight samples (2.7%) (with cycle thresh- 
old [Cj.] values ranging from 25.4 to 36.9) were found to contain 
influenza virus RNA following a real-time reverse transcription 
(RT) -PGR assay for the influenza matrix gene. AIV was confirmed 
by sequence and BLAST analysis of at least one gene sequence 
(hemagglutinin [HA], nonstructural [NS], or matrix [M]) (Ta- 




FIG 1 Map of sampling sites in Antarctica. Adelie penguins were sampled from two locations in Antarctica: Admiralty Bay, King George Island; and Kopaitik 
Island, Rada Covadonga, Antarctic Peninsula. Inset: An Adelie penguin. The map is courtesy of Marco ViUarroel, reproduced with permission. 
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FIG 2 Phylogenetic tree of Hll (A) and N2 (B) sequences. Sequences from 
the viruses isolated from Adelie penguins in Antarctica are highly diverged 
from other circulating AIVs. The HA phylogenetic tree was constructed using 
maximum likelihood methodology and shows ancestry with Eurasian avian 
Hll viruses, while the NA sequence is highly diverged from all other clades but 
with shared ancestry with North American N2 viruses. Detailed phylogenetic 
trees with virus designations can be found in Fig. S2 in the supplemental 
material. 

ble 1). Six AlV-positive samples were from adult Adelie penguins 
and two from Adelie chicks, with all samples except for one col- 
lected from Rada Covadonga (Table 1 ) . All eight influenza viruses, 
including the single virus from Admiralty Bay, had a high degree 
of homology (>99.7% nucleotide sequence identity) between re- 
spective genes. Following inoculation of these viruses into embry- 
onated hens' eggs, four viruses were successfully cultured (all from 
specimens with a Cj- value of <33 based on the matrix gene real- 
time RT-PCR assay). Full genome analysis of the four isolated 
viruses showed that they were all H11N2 influenza A viruses (Gen- 
Bank sequence accession numbers KJ729348 to KJ729379). These 
four viruses were designated A/ Adelie penguin/Antarctica/178/ 
2013, A/Adelie penguin/ Antarctica/184/2013, A/Adelie penguin/ 
Antarctica/226/2013, A/Adelie penguin/Antarctica/270/2013 (ab- 
breviated as A/peng/Ant/178/13, A/peng/Ant/184/13, A/peng/ 
Ant/226/13, and A/peng/ Ant/270/13, respectively). 

Evolutionary relationship and divergence time estimates. 
Full genome sequences of the four isolated viruses were compared 
with all available animal and human influenza virus sequences in 
public databases to determine the phylogenetic relationships and 
the evolutionary origin of each gene segment (Fig. 2 and 3). De- 
tailed phylogenetic trees with strain designations are shown in 
Fig. S2 and S3 in the supplemental material. Unexpectedly, each 
AIV gene segment from the Adelie penguin viruses was highly 
different from that of contemporary AIVs circulating either in the 
Northern or Southern hemispheres (see Fig. S2 and S3 and Ta- 
ble SI in the supplemental material). Four of the gene segments 
(HA, PB2, nucleoprotein [NP] , and NS) were most closely related 
to North American avian lineage viruses from the 1960s to 1980s 
(90 to 96% identity) (Fig. 2 and 3; see also Table SI), while the N2 
neuraminidase gene formed a basal position to the North Ameri- 
can avian influenza lineage, indicating an early divergence from 



the known AIV gene pool (Fig. 2; see also Fig. S2). The MP gene of 
the penguin virus formed a close relationship with North Ameri- 
can viruses from 1974 to 2001 but appears to be the closest ances- 
tor to a large number of South American AIVs from Chile, Argen- 
tina, and Brazil (Fig. 3; see also Fig. S3E). The PA segments of the 
penguin AIVs also have a close ancestral relationship with South 
American viruses, which together appear on their own phyloge- 
netic clade sharing ancestry with equine influenza H3N8 viruses 
from the 1960s (Fig. 3; see also Fig. S3C). Interestingly, the PBl 
segments of the penguin AIVs are also most closely related to 
equine H3N8 viruses, with significant genetic distance from North 
American AIVs (Fig. 3; see also Fig. S3B). 

Using a molecular clock analysis that incorporates the evolu- 
tionary rate of each AIV gene segment, we estimated the time of 
the most recent common ancestor (TMRCA) of the penguin virus 
and its most closely related sequence. For the HA, NA, PB2, NP, 
PA, and NS genes, the TMRCA was between 1957 and 1964, 
whereas for the MP and PBl genes, the estimated TMRCA was 
even earlier at 1933 and 1936, respectively (Table 2). Therefore, 
there has been an estimated 49 to 80 years of unsampled diversity 
and evolution leading to these previously unknown influenza vi- 
ruses (Table 2). 

Molecular markers of host specificity and pathogenicity. The 

HA of the Adelie penguin viruses did not contain the multibasic 
cleavage site observed in highly pathogenic H5 or H7 AIVs, nor 
did the genome contain other major mammalian adaptations, 
such as the PB2 E627K substitution, which is known to alter host 
range (17) and virulence (18). Infection of embryonated eggs was 
not pathogenic to the embryo. The receptor binding pocket resi- 
dues have only one human-like amino acid at position A137 (H3 
numbering), with the rest being avian-like with E190, G225, 
Q226, and G228, suggesting an Q!2,3-linked sialic acid receptor 
binding preference (see Fig. S4 in the supplemental material). 
Compared with the distantly related HI 1N9 virus isolated from a 
wild bird in Australia in 2004, the A/peng/Ant/178/13 HA se- 
quence contained amino acid substitutions at I155T (at the base of 
the receptor binding site and within contact range to the sialic 
acid) and A186S and T187S (at the entry of the host receptor 
pocket), which may alter preferences for different variants of 2,3- 
linked glycans (Fig. 4). Furthermore, A/peng/Ant/178/13 had a 
deletion of the A 134 residue compared to other wild bird Hll 
viruses and an insertion of G133 compared to H13 and H16 vi- 
ruses with known crystal structures (see Fig. S4), which could alter 
receptor binding slightly due to relative shifts of the adjacent 160 
loop. Even though their PBl and PA genes were phylogeneticaUy 
related to those of the equine viruses, the Adelie penguin viruses 
contained typical avian-like amino acid signatures (N375 [PBl], 
D55 and VlOO [PA] ), which differ from the equine-like signatures 
that presumably arose following adaptation to viral replication in 
horses (19). The viruses contained none of the mutations known 
to confer resistance to either the adamantane or neuraminidase 
inhibitor classes of influenza antiviral drugs. 

Identification of influenza antibodies in penguins. Sera from 
270 of the 301 penguins, including the eight penguins that were 
AIV positive, were tested by enzyme-linked immunosorbent assay 
(ELISA) for cross-reactive NP antibodies. As the NP is highly con- 
served across all influenza A viruses, the assay is designed to detect 
any previous influenza A infection, regardless of subtype. In com- 
parison, a hemagglutination inhibition (HI) assay, which detects 
antibodies to the highly variable HA protein, was used to detect 
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Antarctica 

South American avian 
otiier avian 
equine H3N8 
other mammalian 



FIG 3 Phylogenetic trees of the internal gene sequences of avian, equine, and mammalian influenza A viruses. Major avian influenza clades are broadly classified 
into Western, Eastern, or global based on the geographic location of sample collection. Phylogenetic trees were constructed using maximum likelihood 
methodology. AIV sequences from the penguins are marked as "Antarctica," with branches colored red. Blue branches indicate South American avian influenza 
sequences, green branches indicate avian influenza viruses from other regions, and orange branches indicate equine influenza viruses. Detailed phylogenetic trees 
with virus designations can be found in Fig. S3 in the supplemental material. 
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TABLE 2 Summaiy of evolutionary relationships of each gene segment of influenza viruses isolated from Adelie penguins in Antarctica during 
2013 





% identity"^ 


Lineage^ 


Mean time of unsampled diversity in yrs (BCI)"" 


HA gene 


90 


North American avian 


50 (43-58) 


NA gene 


90 


North American avian/ new lineage 


54 (50-60) 


PB2 gene 


96 


North American avian 


50 (41-58) 


PBl gene 


89 


Equine H3N8 


77 (61-100) 


PA gene 


91 


Equine H3N8, South American avian 


56 (50-69) 


NP gene 


96 


North American avian 


49 (41-58) 


MP gene 


97 


North American avian, South American avian 


80 (56-112) 


NS gene 


95 


North American avian 


56 (45-69) 



" Percent identity was estimated using BLAST in NCBI influenza virus resources. 

Virus lineages were determined using a maximum likelihood phylogeny generated for each gene segment of all influenza virus sequences that were >96% similar. Phyiogenies 
with virus names and accession numbers are presented in Fig. SIA to H in the supplemental material. 

Time tothemostrecentcommonancestor(in years ) of influenza viruses isolated from Adelie penguins in Antarctica during 20 1 3 and the most closely related virus sequence; 
represents the mean time of unsampled diversity of this virus (orange bars in Fig. S2 and S3 in the supplemental material}. Bayesian confidence intervals (BCI) for the time of 
unsampled diversity are shown in parentheses. 



previous infection with either the penguin H11N2 viruses or a 
distantly related H11N9 virus isolated from a wild bird in Austra- 
ha in 2004. A total of 15.9% of penguins (43/270) were shown to 
have influenza A NP antibodies by ELISA. All NP antibody- 
positive penguins were adults, with a greater proportion from 




FIG 4 HA structural homology modeling. The structures of penguin and a 
related sandpiper HI 1 virus were modeled based on the HA crystal structures 
of H13 (PDB ID number 4kps), H16 (PDB ID number 4f23), and H2 (PDB ID 
number 2wr3) in complex with an a2,3Gal-linked sialic acid analogue (avian- 
like receptor) . The receptor binding pocket residues for both H 1 1 s are typically 
avian, with E190, G225, Q226, and G228 (H3 numbering), suggesting a2,3- 
linked SA receptor (avian-like) binding (purple and green balls). Comparison 
of the receptor binding pocket in the A/peng/Ant/178/13 with that in A/sharp 
tailed sandpiper/ Australia/10/2004 reveals differences at 1155, A186, and T187 
(red or blue sticks for penguin or sandpiper Hll, respectively), which may 
influence preferences for different variants of 2,3-linked glycans. Several major 
antigenic sites also differ among the HI Is; for example, epitope Sb (at the top 
of the 190 helix) and upper parts of Sa (the area around the 160 loop). Another 
unique structural feature distinguishing the penguin Hll from the nearest 
known Hll sequences and the H13 and H16 structures are different indel 
patterns around G133 (yellow circle). 



Rada Covadonga than from Admiralty Bay (27.3% compared to 
16.7%) (Table 1). Only 2/270 (0.7%) penguins had positive HI 
titers against the penguin H11N2 virus (Table 1). One of the 
HI 1N2 antibody-positive serum samples (with a titer of 40) was 
from an H11N2 PGR- and isolation-positive penguin that yielded 
A/peng/Ant/226/13, while the other positive serum sample (with a 
titer of 80 to 160) was from a penguin that was RT-PCR negative 
for influenza virus. Neither these sera nor the remaining 268 sera 
inhibited the unrelated H11N9 wild bird virus from Australia. 
Overall, 5 of the 8 influenza virus-positive penguins had high NP 
antibody titers, but only one had specific Hll antibodies. This 
suggests that the detected NP antibodies may have been induced 
by prior infection(s), as the remaining three influenza virus- 
positive penguins (which included samples with the highest RNA 
copy number and two that were culture positive) were NP anti- 
body negative. 

Infectivity in ferrets and receptor binding. To assess the po- 
tential risk of A/peng/ Ant/270/13 infection in mammals, we ad- 
ministered the virus to ferrets via intranasal (i.n.) or intratracheal 
(i.t.) routes and compared infectivity with an Australian wild bird 
HI 1N9 virus and a contemporary human H3N2 influenza virus. 
Intranasal administration of A/peng/ Ant/270/1 3 resulted in little 
or no viral replication, with detection of virus only on day 2 
postinfection (at a low titer of 1 X 10' 50% egg infectious doses 
[EID5o]/ml) (see Fig. S5 in the supplemental material). The ferrets 
showed no clinical symptoms or significant weight loss and had no 
detectable Hll serum antibodies by HI assay on day 14 or at 
21 days following an intramuscular boost on day 14. In compari- 
son, i.n. delivery of an Australian wQd bird Hll virus and the 
human H3N2 virus yielded productive infections (shedding virus 
from days 2 to 5, with peak viral titers of 1 X lO''-^ and 1 X 10^-^ 
EID5o/ml, respectively) (see Fig. S5) and homologous HI antibody 
titers of 40 and 2,560, respectively, on day 14 postinfection. Fol- 
lowing i.t. delivery of A/peng/Ant/270/13, a low virus titer (1 X 
10' EIDjo/ml) was detected in one lung lobe but not in the other 
lobe at day 4, but again no HI 1 antibodies were detected on day 
14. No viral titer was detected in either lung lobe following i.t. 
administration of ferrets with either the Australian wild bird Hll 
virus or a human H3N2 influenza virus (see Fig. S5). Using turkey 
erythrocytes enzymaticaUy modified to contain either terminal 
SAQ!2,3Gal (avian-like receptors) or SAa2,6Gal (human-like re- 
ceptors) linkages on cell surface oligosaccharides, we further dem- 
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onstrated that the A/peng/ Ant/2 70/13 primarily bound avian-like 
receptors, with only minor binding to a2,6-linked SA receptors 
(see Table S2). We also showed that, compared to the human 
H3N2 virus, the viral infectivity of A/peng/ Ant/270/ 1 3 in vitro was 
considerably reduced at 33°C compared to that at 37°C, indicating 
a lack of adaptation to replication in the cooler mammalian upper 
airways (see Fig. S6). 

DISCUSSION 

Here, we describe the detection of AIVs in Antarctica, providing 
the first insights into their ecology in this isolated continent. Pre- 
vious studies of birds and mammals in Antarctica have failed to 
detect either live AIV or influenza virus RNA. Although AIV an- 
tibodies have been detected in penguins at frequencies ranging 
from 11.8 to 12.5% in Gentoo penguins {Pygoscelis papua) (12, 
15), from 3.7 to 58.4% in Adelie penguins (12, 14, 16), and 10.8% 
in chinstrap penguins (Pygoscelis antarcticus) (12), several other 
studies have failed to detect any influenza antibodies in these or 
other penguin species in Antarctica (14, 16, 20, 21), suggesting 
that influenza infections may be absent from penguin colonies in 
Antarctica for long periods of time. The detection of HI 1N2 virus 
in 2.7% of Adelie penguins sampled in this study, along with the 
substantially lower prevalence of penguins with detectable Hll 
antibody titers (0.7%), suggests either that we sampled during the 
acute phase of infection before antibodies to HI 1 had developed 
or that the H11N2 infections caused only a weak antibody re- 
sponse. The intensity and longevity of antibody responses in pen- 
guins infected with AIVs are currently unknown, highlighting the 
need for longitudinal repeat sampling studies of the same penguin 
colonies with tagged individuals to enable interpretation of serol- 
ogy data in these species. Although Hll AIVs have never been 
reported in wild birds from South America, they have been de- 
tected in migratory shorebirds in North America (where they were 
the second most frequently detected HA subtype) (22), as well as 
in Australia (23) and New Zealand (24). Additionally, Hll AIVs 
have been commonly detected in ducks in many other regions of 
the world (25, 26). 

It is important to understand the frequency and magnitude of 
inter- and intracontinental movement of Antarctic AIVs. Long- 
distance migrants such as skuas and giant petrels that traverse 
large distances both within Antarctica and between continents 
have been shown to have a high frequency of AIV antibodies 
(skuas, 1.0 to 11.1% [12-14]; giant petrels, 51.5 to 58.8% [12]). 
These birds have very close interactions with penguins, often nest- 
ing alongside their colonies and attacking penguin chicks, provid- 
ing a likely route for AIV transmission. The yellow-bUled pintail 
{Anas georgica) is a species of duck that, while most commonly 
observed in South America, has occasionally been recorded on the 
Antarctic peninsula, making it the most southerly recorded water- 
fowl species (27). Given the high rates of AIV carriage in waterfowl 
compared to those in other orders of birds, it is possible that this 
species plays a role in AIV ecology in Patagonia, the South Shet- 
land Islands, and the Antarctic Peninsula. The return of migratory 
birds to Antarctica each year between October and December 
(start of the austral summer period) coincides with penguins re- 
turning to colonies to nest and rear chicks. Because of the high 
population density of birds, summer is the most likely period for 
AIV transmission and infection in Antarctica. The large amount 
of penguin feces in colonies during summer, which in some cases 
is so significant it can be observed on satellite images (28), pre- 



sumably facilitates AIV transmission by the fecal-oral route (29). 
Although the relative susceptibility of different penguin species to 
AIVs is unknown, we did take swabs from a small number of 
gentoo (« = 74) and chinstrap (n = 18) penguins in Rada Cova- 
donga, but none were AIV positive. The different species typically 
form exclusive separate colonies, which may explain the lack of 
cross-infection between the three penguin species at the time of 
sampling in Rada Covadonga. 

The detection of this highly divergent H 1 1N2 AIV suggests that 
its existence in Antarctica is not transient. Although two gene 
segments share common ancestry with South American AIVs, the 
other remaining six genes shared ancestry with either North 
American avian or equine viruses. Antarctica therefore maybe the 
most extreme example of an AIV "evolutionary sink," where vi- 
ruses are seeded into the continent on an infrequent basis and then 
become established, persisting as a local reservoir and evolving 
independently from AIVs in other areas of the globe due to a lack 
of interconnectedness. Future studies to understand which species 
or groups of animals are responsible for maintaining AIVs in Ant- 
arctica are needed. Based on all previous studies in Antarctica, the 
overall AIV antibody prevalence in penguins (3.3%; 82/2,500) (12, 
14-16, 20, 21, 30, 31) is considerably lower than it is in skuas and 
giant petrels (13.3%; 47/353) (12-14), suggesting that the latter 
species may be more likely to maintain and transmit AIVs in Ant- 
arctica. An alternative hypothesis is that AIVs circulate between 
birds only in summer and then become cryopreserved in ice over 
winter (32), after which birds are again infected once the ice melts 
in the following summer. In addition, a number of Antarctic ma- 
rine mammals, such as seals and whales, may also have the poten- 
tial to become infected and spread influenza viruses given the 
previous reports of their susceptibility to influenza infection (33, 
34). Although three reported studies of seals (n = 306) in Antarc- 
tica did not detect influenza viruses or antibodies (14, 35, 36), 
these animals appear to be susceptible to both avian and human 
influenza viruses, as their respiratory tracts contain both q:2,3- 
and a2,6-linked sialic acid receptors (33). Seals have been infected 
with both AIV (33) and human influenza viruses (37), as well as 
infecting humans with an H7N7 virus following close experimen- 
tal contact (38, 39). Whether these animals could act as a "mixing 
vessel" for avian or human influenza viruses in a natural environ- 
ment in the same manner as pigs is unknown, although the bidi- 
rectional transmission of influenza between seals and humans 
may be a potential risk for both local wildlife and public health 
with increasing human activity in Antarctica. The HI 1N2 Antarc- 
tica AIV A/peng/Ant/270/13 appeared to be less infectious in fer- 
rets than has been observed for other AIVs from wild birds. AIVs 
such as H6N1, H1N9, H9N2, H6N5, and H11N9 (used in this 
study) have readily infected ferrets and in many cases been able to 
successfully transmit between cohoused naive ferrets (40-42). 
Further work is needed to understand the key residues that may be 
responsible for these differences in replication. 

While it has previously been postulated that AIV infection 
caused mortality of Adelie penguin chicks in Antarctica (16), we 
did not notice any clinical symptoms or below-average weight in 
the HllN2-infected penguins in this study. Given that highly 
pathogenic AIVs have been detected in both South America (43) 
and Australia (44), it is not inconceivable that such viruses could 
be transferred to the Antarctic continent by migratory birds, po- 
tentially resulting in catastrophic mass mortality, such as that seen 
with other pathogens such as avian cholera (Pasteurella multocida) 
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(45, 46). This is the first insight into AIV ecology in Antarctica and 
provides further evidence, in addition to some South American 
studies (8, 47, 48), of an avian ancestry to H3N8 equine influenza 
viruses that emerged and spread globally among horses since 1963 
(49). As such, this study and the recent detection of novel influ- 
enza A viruses in fruit bats from Central and South America (50, 
51) fill important gaps in our understanding of the global move- 
ment of AIVs, the reservoir species that may maintain them, and 
their potential impact on animal and human health. 

MATERIALS AND METHODS 

Sampling locations and times. Adelie penguins were sampled from two 
locations in Antarctica (Fig. 1): Admiralty Bay, King George Island 
(62°9'35"S, 58°28'17"W), 1 km east of Arctowski Polish Base, between 
14 January and 31 January 2013; and Kopaitik Island, Rada Covadonga, 
1 km west of Base General Bernardo O'Higgins Riquelme, Antarctic Pen- 
insula (63°19'15"S, 57°53'55"W), between 1 Februaiy and 15 February 
2013. 

Sample collection from penguins. Separate cloacal and tracheal swabs 
were taken from each Adelie penguin and then placed in the same tube 
containing viral transport medium (brain heart infusion [BHI] broth- 
based medium [Oxoid; catalog number CM1135] with 0.3 mg/ml peni- 
cillin, 5 mg/ml streptomycin, 0.1 mg/ml gentamicin, and 2.5 p.g/ml am- 
photericin B). Swab specimens were kept on ice for up to 4 h before being 
frozen at — 80°C. Blood was collected from penguins either from wing or 
foot veins and put into serum clot- activating tubes, centrifuged, and se- 
rum stored at either — 80°C or — 20°C. Approvals to conduct sampling 
fi-om penguins in Antarctica were provided by the Universidad de Con- 
cepcion, Facultad de Ciencias Veterinarias, Chilian, Chile (application 
number CE-3-2010), and Instituto Antartico Chileno, Chile (application 
number 03/2013). 

RT-PCR and sequence analysis. RNA was isolated from cloacal/tra- 
cheal swabs using the QIAGEN QIAextractor system. AU samples were 
tested by a one-step RT-PCR assay targeting the influenza matrix genome 
segment (52), using a SensiFast Probe Lo-ROX one-step kit (Bioline) on 
an Applied Biosystems 7500 Fast real-time PGR system (Life Technolo- 
gies). Samples with a cycle threshold (Cj.) value of <42 cycles or cultured 
influenza viruses were analyzed further by conventional RT-PCR using 
universal primers targeting each of the eight influenza segments (53). 
Amplified RT-PCR products were sequenced using standard techniques 
on an Applied Biosystems 3500XL genetic analyzer. 

Virus isolation. Virus isolation was attempted from RT-PCR-positive 
samples by inoculating the swab specimen (diluted 1:1 with phosphate- 
buffered saline [PBS] containing 1% neomycin-polymyxin solution [bi- 
oCSL]) into the allantoic cavity of 11 -day-old embryonated hens' eggs. 
Allantoic fluid was harvested after 3 days, and influenza virus was detected 
by hemagglutination with turkey erythrocytes. The EIDgg/ml of each iso- 
late was determined by performing logi^ dilutions in PBS containing 1% 
neomycin-polymyxin solution and infecting the allantoic cavity of 1 1 - 
day-old embryonated hens' eggs. Calculation of infectious titers was per- 
formed according to the Reed and Muench method (54). 

Serological analysis. To test serum samples for broadly reactive influ- 
enza NP antibodies, a competitive ELISA with plates coated with the NP 
from A/ turkey/Ontario/62 13/66 (H5N1) was used as described by Selleck 
and Kirkland (55). Serum was considered antibody positive if the sample 
inhibited >60% of the monoclonal antibody (MAb) binding. To test se- 
rum samples for specific HA antibodies, a hemagglutination inhibition 
(HI) assay using turkey erythrocytes was used following previously de- 
scribed methods (56). 

Dataset preparation and preliminary phylogenetic analysis. To elu- 
cidate the evolutionary relationships of viruses isolated in this study, we 
first carried out a provisional phylogenetic analysis for the hemagglutinin 
and neuraminidase genes, including all available Hll and N2 subtype 
influenza gene sequences, and for the internal gene datasets, including 
viruses of all subtypes, using the maximum likelihood method with a 



general time reversible (GTR) nucleotide substitution model in FastTree 2 
(57). The purpose of these large-scale phylogenetic analyses was to iden- 
tify broad relationships between the new sequences and established influ- 
enza A viruses from multiple hosts. These analyses indicated that the 
Adelie penguin viruses were distantly related to previously detected avian 
influenza virus and their derived lineages; hence, for all subsequent anal- 
yses, the dataset sizes were reduced by identifying representative se- 
quences at 96% identity using the software CD-HIT (58). Due to the 
geographical proximity of Antarctica to South America, sequences of all 
publicly available avian isolates from South America were included in our 
final datasets. 

Maximum likelihood and Bayesian relaxed-molecular clock analy- 
ses. For the final datasets, maximum likelihood phylogenies were gener- 
ated using the GTR model with gamma rate heterogeneity, using RaXML 
(59). Statistical supports were estimated using 500 maximum likelihood 
bootstrap replicates. To estimate divergence times, we applied the log- 
normal relaxed-clock Bayesian Markov chain Monte Carlo method with 
an SRD06 codon substitution model (60), as implemented in BEAST ver- 
sion 1.7.5 (61). For each dataset, multiple analyses were conducted for 50 
million generations, sampling every 5,000 generations, resulting in 10,000 
trees. After convergence was ensured, the final 90% of trees were summa- 
rized using TreeAnnotator (61) and produced using Fig tree (http://tree 
.bio.ed.ac.uk/software/figtree/). 

HA structural homology modeling. Due to the lack of a published 
Hll crystal structure, numerous models were made for the Hll viruses 
based on the crystal structure of the most closely related HAs available in 
the RCSB protein databank (http://www.rcsb.org). The Hll structures 
were modeled based on the HA crystal structures of H 13 (PDB ID number 
4kps), H16 (PDB ID number 4f23), and H2 (PDB ID number 2wr3) in 
complex with an a2,3Gal-linked sialic acid analogue (avian-like receptor) 
using MODELLER (62) and YASARA (63). 

Ferret infectivity. Ferrets 6 to 18 months old were administered 1 X 
10^ EID5o/ml egg-propagated virus via the i.n. or i.t. routes. Of two ferrets 
administered virus via the i.n. route, one was nasal washed on days 2,3,4, 
and 5 before being bled and euthanized on day 14, and the other was bled 
and given a subsequent intramuscular boost with inactivated virus on day 
14 and then bled again and euthanized on day 21. Of two ferrets admin- 
istered virus via the i.t. route, one was nasal washed on days 2, 3, 4, and 5 
before being bled and euthanized on day 14, and the other ferret was 
culled on day 4 to analyze viral titers in lungs. Weights, temperatures, and 
clinical signs were monitored throughout the experiment. 

Receptor binding assay. The receptor specificity of viruses was as- 
sessed using methods described by Rogers and Paulson (64). Briefly, vi- 
ruses were diluted to a hemagglutination titer of 64 hemagglutinating 
units (HAU)/25 jA in PBS. A 10% suspension of turkey erythrocytes was 
treated for 1 h at 37°C with 500 mU/ml Vibrio cholerae enzyme (Sigma). 
Erythrocytes were washed three times with PBS and resuspended to a 20% 
suspension. The erythrocytes were then mixed with an equal volume of 
recombinant beta-galactoside-alpha-2,3-sialtransferase (10 mU//j.l) or 
recombinant beta-galactoside-alpha-2,6-sialtransferase (10 mU/p,l) (Ja- 
pan Tobacco, Inc.) and incubated for 3 to 4 h at 37°C. Erythrocytes were 
washed a further three times with PBS and resuspended to 1%. Hemag- 
glutination assays were conducted using untreated erythrocytes, sialidase- 
treated erythrocytes (Asialo), and sialidase-treated erythrocytes resialy- 
lated to contain either terminal SAa2,3Gal or SAQ!2,6Gal linkages on cell 
surface oligosaccharides. 

Cell culture infectivity assay at 33°C and 37°C. Madin-Darby canine 
kidney (MDCK) cells (American Type Culture Collection; CCL-34) were 
grown to near confluence in 96-well plates (5% CO2, 37°C). LogjQ virus 
dilutions were made in Dulbecco's modified Eagle's medium (SAFC Bio- 
sciences) containing sodium bicarbonate (3%) with the addition of 2 mM 
glutamine (SAFC Biosciences), 1% nonessential amino acids (Sigma), 
0.05% NaHC03 (Sigma), 0.02 M HEPES (Sigma), 4% penicillin and 
streptomycin (Sigma), and 4 fJ.g/ml trypsin (Sigma). Cells were washed 
twice with PBS before infection with four replicates of 10-fold dilutions of 
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each virus. After incubation for 5 days at 33°C or 37°C in 5% COj, each 
well was scored for virus growth by cytopathic effect and hemagglutina- 
tion of turkey erythrocytes, and the dose required to infect 50% of wells 
(50% tissue culture infectious dose [TCIDgg]) was determined by the 
Reed and Muench method (54). 

Nucleotide sequence accession numbers. FuU genome sequences 
from the four H11N2 influenza viruses isolated in eggs have been depos- 
ited in GenBank under accession numbers KJ729348 to KJ729379. 
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